PCR candidate region mismatch scanning: adaptation to quantitative, high-throughput genotyping

Linkage and association analyses were performed to identify loci affecting disease susceptibility by scoring previously characterized sequence variations such as microsatellites and single nucleotide polymorphisms. Lack of markers in regions of interest, as well as difficulty in adapting various methods to high-throughput settings, often limits the effectiveness of the analyses. We have adapted the Escherichia coli mismatch detection system, employing the factors MutS, MutL and MutH, for use in PCR-based, automated, high-throughput genotyping and mutation detection of genomic DNA. Optimal sensitivity and signal-to-noise ratios were obtained in a straightforward fashion because the detection reaction proved to be principally dependent upon monovalent cation concentration and MutL concentration. Quantitative relationships of the optimal values of these parameters with length of the DNA test fragment were demonstrated, in support of the translocation model for the mechanism of action of these enzymes, rather than the molecular switch model. Thus, rapid, sequence-independent optimization was possible for each new genomic target region. Other factors potentially limiting the flexibility of mismatch scanning, such as positioning of dam recognition sites within the target fragment, have also been investigated. We developed several strategies, which can be easily adapted to automation, for limiting the analysis to intersample heteroduplexes. Thus, the principal barriers to the use of this methodology, which we have designated PCR candidate region mismatch scanning, in cost-effective, high-throughput settings have been removed.

A quantitative, high-throughput screen for protein stability

In proteomic research, it is often necessary to screen a large number of polypeptides for the presence of stable structure. Described here is a technique (referred to as SUPREX, stability of unpurified proteins from rates of H/D exchange) for measuring the stability of proteins in a rapid, high-throughput fashion. The method uses hydrogen exchange to estimate the stability of microgram quantities of unpurified protein extracts by using matrix-assisted laser desorption/ionization MS. The stabilities of maltose binding protein and monomeric �� repressor variants determined by SUPREX agree well with stability data obtained from conventional CD denaturation of purified protein. The method also can detect the change in stability caused by the binding of maltose to maltose binding protein. The results demonstrate the precision of the method over a wide range of stabilities.

An automated multiplex oligonucleotide synthesizer: development of high-throughput, low-cost DNA synthesis.

An automated oligonucleotide synthesizer has been developed that can simultaneously and rapidly synthesize up to 96 different oligonucleotides in a 96-well microtiter format using phosphoramidite synthesis chemistry. A modified 96-well plate is positioned under reagent valve banks, and appropriate reagents are delivered into individual wells containing the growing oligonucleotide chain, which is bound to a solid support. Each well has a filter bottom that enables the removal of spent reagents while retaining the solid support matrix. A seal design is employed to control synthesis environment and the entire instrument is automated via computer control. Synthesis cycle times for 96 couplings are < 11 min, allowing a plate of 96 20-mers to be synthesized in < 5 hr. Oligonucleotide synthesis quality is comparable to commercial machines, with average coupling efficiencies routinely > 98% across the entire 96-well plate. No significant well-to-well variations in synthesis quality have been observed in > 6000 oligonucleotides synthesized to date. The reduced reagent usage and increased capacity allow the overall synthesis cost to drop by at least a factor of 10. With the development of this instrument, it is now practical and cost-effective to synthesize thousands to tens of thousands of oligonucleotides.

SpliceDB: database of canonical and non-canonical mammalian splice sites

A database (SpliceDB) of known mammalian splice site sequences has been developed. We extracted 43 337 splice pairs from mammalian divisions of the gene-centered Infogene database, including sites from incomplete or alternatively spliced genes. Known EST sequences supported 22��815 of them. After discarding sequences with putative errors and ambiguous location of splice junctions the verified dataset includes 22 489 entries. Of these, 98.71% contain canonical GT���AG junctions (22 199 entries) and 0.56% have non-canonical GC���AG splice site pairs. The remainder (0.73%) occurs in a lot of small groups (with a maximum size of 0.05%). We especially studied non-canonical splice sites, which comprise 3.73% of GenBank annotated splice pairs. EST alignments allowed us to verify only the exonic part of splice sites. To check the conservative dinucleotides we compared sequences of human non-canonical splice sites with sequences from the high throughput genome sequencing project (HTG). Out of 171 human non-canonical and EST-supported splice pairs, 156 (91.23%) had a clear match in the human HTG. They can be classified after sequence analysis as: 79 GC���AG pairs (of which one was an error that corrected to GC���AG), 61 errors corrected to GT���AG canonical pairs, six AT���AC pairs (of which two were errors corrected to AT���AC), one case was produced from a non-existent intron, seven cases were found in HTG that were deposited to GenBank and finally there were only two other cases left of supported non-canonical splice pairs. The information about verified splice site sequences for canonical and non-canonical sites is presented in SpliceDB with the supporting evidence. We also built weight matrices for the major splice groups, which can be incorporated into gene prediction programs. SpliceDB is available at the computational genomic Web server of the Sanger Centre: http://genomic.sanger.ac.uk/spldb/SpliceDB.html and at http://www.softberry.com/spldb/SpliceDB.html.

Alkaline-mediated differential interaction (AMDI): A simple automatable single-nucleotide polymorphism assay

The key requirements for high-throughput single-nucleotide polymorphism (SNP) typing of DNA samples in large-scale disease case-control studies are automatability, simplicity, and robustness, coupled with minimal cost. In this paper we describe a fluorescence technique for the detection of SNPs that have been amplified by using the amplification refractory mutation system (ARMS)-PCR procedure. Its performance was evaluated using 32 sequence-specific primer mixes to assign the HLA-DRB alleles to 80 lymphoblastoid cell line DNAs chosen from our database for their diversity. All had been typed previously by alternative methods, either direct sequencing or gel electrophoresis. We believe the detection system that we call AMDI (alkaline-mediated differential interaction) satisfies the above criteria and is suitable for general high-throughput SNP typing.

Direct genetic analysis by matrix-assisted laser desorption/ionization mass spectrometry

An approach to analyzing single-nucleotide polymorphisms (SNPs) found in the human genome has been developed that couples a recently developed invasive cleavage assay for nucleic acids with detection by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). The invasive cleavage assay is a signal amplification method that enables the analysis of SNPs by MALDI-TOF MS directly from human genomic DNA without the need for initial target amplification by PCR. The results presented here show the successful genotyping by this approach of twelve SNPs located randomly throughout the human genome. Conventional Sanger sequencing of these SNP positions confirmed the accuracy of the MALDI-TOF MS analysis results. The ability to unambiguously detect both homozygous and heterozygous genotypes is clearly demonstrated. The elimination of the need for target amplification by PCR, combined with the inherently rapid and accurate nature of detection by MALDI-TOF MS, gives this approach unique and significant advantages in the high-throughput genotyping of large numbers of SNPs, useful for locating, identifying, and characterizing the function of specific genes.

Quantitative high performance liquid chromatography/tandem mass spectrometric analysis of the four classes of F2-isoprostanes in human urine

Isoprostanes (iPs) are free radical catalyzed prostaglandin isomers. Analysis of individual isomers of PGF2�����F2-iPs���in urine has reflected lipid peroxidation in humans. However, up to 64 F2-iPs may be formed, and it is unknown whether coordinate generation, disposition, and excretion of F2-iPs occurs in humans. To address this issue, we developed methods to measure individual members of the four structural classes of F2-iPs, using liquid chromatography/tandem mass spectrometry (LC/MS/MS), in which sample preparation is minimized. Authentic standards of F2-iPs of classes III, IV, V, and VI were used to identify class-specific ions for multiple reaction monitoring. Using iPF2��-VI as a model compound, we demonstrated the reproducibility of the assay in human urine. Urinary levels of all F2-iPs measured were elevated in patients with familial hypercholesterolemia. However, only three of eight F2-iPs were elevated in patients with congestive heart failure, compared with controls. Paired analyses by GC/MS and LC/MS/MS of iPF2��-VI in hypercholesterolemia and of 8,12-iso-iPF2��-VI in congestive heart failure were highly correlated. This approach will permit high throughput analysis of multiple iPs in human disease.

Serial microanalysis of renal transcriptomes

Large-scale gene expression studies can now be routinely performed on macroamounts of cells, but it is unclear to which extent current methods are valuable for analyzing complex tissues. In the present study, we used the method of serial analysis of gene expression (SAGE) for quantitative mRNA profiling in the mouse kidney. We first performed SAGE at the whole-kidney level by sequencing 12,000 mRNA tags. Most abundant tags corresponded to transcripts widely distributed or enriched in the predominant kidney epithelial cells (proximal tubular cells), whereas transcripts specific for minor cell types were barely evidenced. To better explore such cells, we set up a SAGE adaptation for downsized extracts, enabling a 1,000-fold reduction of the amount of starting material. The potential of this approach was evaluated by studying gene expression in microdissected kidney tubules (50,000 cells). Specific gene expression profiles were obtained, and known markers (e.g., uromodulin in the thick ascending limb of Henle's loop and aquaporin-2 in the collecting duct) were found appropriately enriched. In addition, several enriched tags had no databank match, suggesting that they correspond to unknown or poorly characterized transcripts with specific tissue distribution. It is concluded that SAGE adaptation for downsized extracts makes possible large-scale quantitative gene expression measurements in small biological samples and will help to study the tissue expression and function of genes not evidenced with other high-throughput methods.

A quantitative method for evaluating the stabilities of nucleic acids

We report a general method for screening, in solution, the impact of deviations from canonical Watson-Crick composition on the thermodynamic stability of nucleic acid duplexes. We demonstrate how fluorescence resonance energy transfer (FRET) can be used to detect directly free energy differences between an initially formed ���reference��� duplex (usually a Watson-Crick duplex) and a related ���test��� duplex containing a lesion/alteration of interest (e.g., a mismatch, a modified, a deleted, or a bulged base, etc.). In one application, one titrates into a solution containing a fluorescently labeled, FRET-active, reference duplex, an unlabeled, single-stranded nucleic acid (test strand), which may or may not compete successfully to form a new duplex. When a new duplex forms by strand displacement, it will not exhibit FRET. The resultant titration curve (normalized fluorescence intensity vs. logarithm of test strand concentration) yields a value for the difference in stability (free energy) between the newly formed, test strand-containing duplex and the initial reference duplex. The use of competitive equilibria in this assay allows the measurement of equilibrium association constants that far exceed the magnitudes accessible by conventional titrimetric techniques. Additionally, because of the sensitivity of fluorescence, the method requires several orders of magnitude less material than most other solution methods. We discuss the advantages of this method for detecting and characterizing any modification that alters duplex stability, including, but not limited to, mutagenic lesions. We underscore the wide range of accessible free energy values that can be defined by this method, the applicability of the method in probing for a myriad of nucleic acid variations, such as single nucleotide polymorphisms, and the potential of the method for high throughput screening.

The Medicago Genome Initiative: a model legume database

The Medicago Genome Initiative (MGI) is a database of EST sequences of the model legume Medicago truncatula. The database is available to the public and has resulted from a collaborative research effort between the Samuel Roberts Noble Foundation and the National Center for Genome Resources to investigate the genome of M.truncatula. MGI is part of the greater integrated Medicago functional genomics program at the Noble Foundation (http://www.noble .org), which is taking a global approach in studying the genetic and biochemical events associated with the growth, development and environmental interactions of this model legume. Our approach will include: large-scale EST sequencing, gene expression profiling, the generation of M.truncatula activation-tagged and promoter trap insertion mutants, high-throughput metabolic profiling, and proteome studies. These multidisciplinary information pools will be interfaced with one another to provide scientists with an integrated, holistic set of tools to address fundamental questions pertaining to legume biology. The public interface to the MGI database can be accessed at http://www.ncgr.org/research/mgi.

trEST, trGEN and Hits: access to databases of predicted protein sequences

High throughput genome (HTG) and expressed sequence tag (EST) sequences are currently the most abundant nucleotide sequence classes in the public database. The large volume, high degree of fragmentation and lack of gene structure annotations prevent efficient and effective searches of HTG and EST data for protein sequence homologies by standard search methods. Here, we briefly describe three newly developed resources that should make discovery of interesting genes in these sequence classes easier in the future, especially to biologists not having access to a powerful local bioinformatics environment. trEST and trGEN are regularly regenerated databases of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Hits is a web-based data retrieval and analysis system providing access to precomputed matches between protein sequences (including sequences from trEST and trGEN) and patterns and profiles from Prosite and Pfam. The three resources can be accessed via the Hits home page (http://hits.isb-sib.ch).

High-sensitivity array analysis of gene expression for the early detection of disseminated breast tumor cells in peripheral blood

Early detection is an effective means of reducing cancer mortality. Here, we describe a highly sensitive high-throughput screen that can identify panels of markers for the early detection of solid tumor cells disseminated in peripheral blood. The method is a two-step combination of differential display and high-sensitivity cDNA arrays. In a primary screen, differential display identified 170 candidate marker genes differentially expressed between breast tumor cells and normal breast epithelial cells. In a secondary screen, high-sensitivity arrays assessed expression levels of these genes in 48 blood samples, 22 from healthy volunteers and 26 from breast cancer patients. Cluster analysis identified a group of 12 genes that were elevated in the blood of cancer patients. Permutation analysis of individual genes defined five core genes (P ��� 0.05, permax test). As a group, the 12 genes generally distinguished accurately between healthy volunteers and patients with breast cancer. Mean expression levels of the 12 genes were elevated in 77% (10 of 13) untreated invasive cancer patients, whereas cluster analysis correctly classified volunteers and patients (P = 0.0022, Fisher's exact test). Quantitative real-time PCR confirmed array results and indicated that the sensitivity of the assay (1:2 �� 108 transcripts) was sufficient to detect disseminated solid tumor cells in blood. Expression-based blood assays developed with the screening approach described here have the potential to detect and classify solid tumor cells originating from virtually any primary site in the body.

Simultaneous assessment of loss of heterozygosity at multiple microsatellite loci using semi-automated fluorescence-based detection: subregional mapping of chromosome 4 in cervical carcinoma.

Detection of loss of heterozygosity (LOH) by comparison of normal and tumor genotypes using PCR-based microsatellite loci provides considerable advantages over traditional Southern blotting-based approaches. However, current methodologies are limited by several factors, including the numbers of loci that can be evaluated for LOH in a single experiment, the discrimination of true alleles versus "stutter bands," and the use of radionucleotides in detecting PCR products. Here we describe methods for high throughput simultaneous assessment of LOH at multiple loci in human tumors; these methods rely on the detection of amplified microsatellite loci by fluorescence-based DNA sequencing technology. Data generated by this approach are processed by several computer software programs that enable the automated linear quantitation and calculation of allelic ratios, allowing rapid ascertainment of LOH. As a test of this approach, genotypes at a series of loci on chromosome 4 were determined for 58 carcinomas of the uterine cervix. The results underscore the efficacy, sensitivity, and remarkable reproducibility of this approach to LOH detection and provide subchromosomal localization of two regions of chromosome 4 commonly altered in cervical tumors.